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Chapter  1 


Introduction 


1.1  Statement  of  the  Problem 

The  research  conducted  here  is  relevant  to  combustion  of  solid  propellants.  The  principal 
objective  of  this  research  is  to  develop  chemical-kinetic  models  which  can  be  employed  to 
model  the  combustion  of  solid  propellants.  The  research  was  performed  in  collaboration  with 
scientists  at  the  U.  S.  Army  Research  Laboratory  at  Aberdeen  Proving  Grounds,  Maryland. 

The  physical  and  chemical  processes  that  govern  the  ignition  and  combustion  of  solid  pro¬ 
pellants  are  complex.  Several  zones  have  been  identified  during  combustion  of  propellants 
(Kubota  1981).  Within  the  propellant,  phase  decomposition  reactions  take  place.  At  the  sur¬ 
face  of  the  regressing  propellant  the  solid  melts  and  then  vaporizes  and  significant  exothermic 
reactions  take  place  in  the  liquid  phase.  In  the  liquid  layer  on  the  surface  of  the  propellant 
bubbles  develop  producing  a  ‘Toam”  reaction  zone.  Above  the  “foam”  zone  many  propellants 
exhibit  a  two-stage  flame  zone.  In  this  two-stage  zone  a  nonluminous  region  separates  the  pri¬ 
mary  reaction  zone  near  the  surface  of  the  propellant  from  the  luminous  secondary  flame  zone. 
The  nonluminous  region  is  commonly  referred  to  as  the  dark  zone.  In  the  primary  reaction  zone 
the  decomposition  of  the  solid  propellant  takes  place  and  moderately  reactive  intermediate 
species  are  formed.  These  intermediate  species  are  convected  away  from  the  surface  and  they 
ignite  after  a  short  period  of  time.  Combustion  then  continues  in  the  secondary  reaction  zone. 
Previous  studies  have  demonstrated  the  importance  of  the  secondary  reaction  zone  in  con¬ 
trolling  the  burning  rate  of  neat  RDX  (Melius  1986,  Melius  1987,  Melius  1988,  Melius  1989). 
Recent  experiments  show  that  chemical  effects  caused  by  the  use  of  propellants  of  differ¬ 
ent  chemical  composition  can  lead  to  major  differences  in  ignition  delays  times  during  typ¬ 
ical  large  calibre  gun  ballistics  cycles  (Kooker,  Chang  &  Howard  1993,  Kooker,  Chang  & 
Howard  1994,  Kooker,  Chang  Sz  Howard  1995,  Kooker  1995).  It  is  believed  that  the  differ¬ 
ences  in  ignition  delay  times  for  different  propellants  arise  from  differences  in  the  chemistry 
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taking  place  in  the  dark  zone  (Kooker  et  al.  1995,  Kooker  1995).  The  focus  of  this  study  is  to 
investigate  the  structure  of  the  secondary  reaction  zone  and  the  structure  of  the  dark  zone. 

The  gas-phase  combustion  of  propellants  takes  place  in  a  number  of  subsystems,  and  the 
reactions  between  CH4-NO2  and  CH2O-NO2  are  such  subsystems.  An  improved  understanding 
of  the  chemistry  taking  place  in  these  systems  can  be  useful  in  modeling  the  combustion  of 
solid  propellants.  A  clear  understanding  of  the  structure  of  hydrocarbon  air  and  alcohol  air 
flames  is  required  to  model  the  structure  of  CH4-NO2  and  CH20-N02  chemical  systems.  The 
research  comprised  several  studies.  They  were  1)  studies  on  the  structures  of  CH4-NO2  and 
CH2O-NO2  flames,  2)  studies  on  the  structures  of  hydrocarbon-air  and  alcohol-air  flames, 
and  3)  studies  on  the  structure  of  the  dark  zone  formed  over  the  burning  surface  of  solid 
propellants.  The  results  of  these  studies  are  summarized  in  the  following  publications  and 
Ph.D  theses. 

1.2  Publications  and  Ph.D  Theses 

1.2.1  Publications 

The  papers  published  and  those  submitted  for  publication  using  support  from  this  research 
program  are 

1.  Bui,  M.,  Seshadri,  K.,  Williams,  F.  A.,  “The  Asymptotic  Structure  of  Premixed  Methane- 
Air  Flames  with  Slow  CO  Oxidation,”  Combust  Flame  89  (1992),  pp  343-362. 

2.  CheUiah,  H.,  Bui-Pham,  M.,  Seshadri,  K.,  and  Law,  C.  K.,  “Numerical  Description  of 
the  Structure  of  Counterflow  Heptane- Air  Flames  Using  Detailed  and  Reduced  Chem¬ 
istry  with  Comparisons  to  Experiments,”  Twenty-Fourth  Symposium  (International)  on 
Combustion^  The  Combustion  Institute,  1992,  pp  851-857. 

3.  Yang,  B.,  and  Seshadri,  K.,  “Asymptotic  Analysis  of  the  Structure  of  Nonpremixed 
Methane  Air  Flames  using  Reduced  Chemistry,”  Combust.  Sci.  and  Tech.  88  (1992), 
pp  115-132. 

4.  Yang,  B.,  Seshadri,  K.,  and  Peters,  N.,  “The  Asymptotic  Structure  of  Premixed  Methanol- 
Air  Flames,”  Combust  Flame  91  (1992),  pp  382-398. 

5.  CheUiah,  H.  K.,  Seshadri,  K.,  and  Law,  C.  K.,  “Reduced  Kinetic  Mechanisms  for  Coun¬ 
terflow  Methane- Air  Diffusion  Flames,”  in  Reduced  Kinetic  Mechanisms  for  Applications 
in  Combustion  Systems  (Eds.  N.  Peters  and  B,  Rogg),  Lecture  Notes  in  Physics^  m  15, 
Springer- Verlag,  Heidelberg,  1993,  pp  224-240. 
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6.  Seshadri,  K.,  and  Williams,  F.  A.,  “Reduced  Chemical  Systems  and  their  Applications 
in  Turbulent  Combustion,”  Chapter  four  in  Turbulent  Reactive  Flows  (Eds.  P,  A.  Libby 
and  F.  A.  Williams),  Academic  Press,  New  York,  1994,  pp  153-210. 

7.  Yang,  B.  and  Seshadri,  K.,  “The  Asymptotic  Structure  of  MethanohAir  Diffusion  Flames,” 
Combust  ScL  and  Tech.  97  (1994),  pp  193-218. 

8.  Trees,  D.,  Brown,  T.  M.,  Seshadri,  K.,  Smooke,  M.  D.,  Balakrishnan,  G.,  Pitz,  R. 
W.,  Giovangigli,  V.,  and  Nandula,  S.  P.,  “The  Structure  of  Nonpremixed  Hydrogen- Air 
Flames,”  Combust.  Sci.  and  Tech.  104  (1995),  pp  427-439. 

9.  Dincic,  N.,  Bui-Pham,  M.  N.,  and  Seshadri,  K.,  “Structure  of  Nonpremixed  CH4/NO2 
Flames,”  submitted  for  publication  in  Combustion  Science  and  Technology,  December 
1995. 

10.  Anderson,  W.  R.,  Dincic,  N.,  and  Seshadri,  K.,  “Studies  of  Reactions  Pertaining  to  and 
Development  of  a  Reduced  Mechanism  for  the  Double  Base  PropeUant  Dark  Zone,”  To 
appear  in  the  proceedings  of  the  31®^  JANNAF  Combustion  Meeting,  Sunnyvale,  CA, 
October,  1994. 

1.2.2  Ph.D  Theses 

The  Ph.D  theses  completed  using  support  from  this  research  program  are 

1.  Mary  N.  Bui-Pham,  “Studies  in  Structures  of  laminar  Hydrocarbon  Flames,”  Ph.D 
Thesis,  University  of  Cahfornia  at  san  Diego,  La  JoUa,  California,  1992. 

2.  Bo  Yang,  “The  Structure  of  Methane,  Methanol  and  Formaldehyde  Flames,”  Ph.D  The¬ 
sis,  University  of  California  at  san  Diego,  La  JoUa,  California,  1993. 

3.  Nenad  Dincic,  “Ignition  and  Extinction  Phenomena  in  Combustion  Systems,”  Ph.D 
Thesis,  University  of  Cahfornia  at  san  Diego,  La  JoUa,  Cahfornia,  1995. 

A  significant  part  of  the  important  results  obtained  in  this  study  are  summarized  in  the 
pubhcations.  Reprints  of  pubhcations  1-8  and  10  are  provided.  Reprints  of  pubhcation  9  wiU 
be  provided  as  soon  as  they  become  available.  A  significant  part  of  the  results  discussed  in  the 
Ph.D  theses  are  also  summarized  in  these  pubhcations.  Chapter  2  summarizes  results  of  our 
studies  on  the  structure  of  nonpremixed  CH4/NO2  flames  which  is  also  included  in  pubhcation 
number  9.  Chapter  3  summarizes  results  of  our  study  on  the  structures  of  CH4/NO2/O2  and 
CH2O/NO2/O2  premixed  flames  which  are  not  available  in  the  pubhcations  and  is  taken  from 
the  thesis  of  Yang  (1993). 
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Chapter  2 


The  Structure  of  Nonpremixed 
CH4/NO2  Flames 


2.1  Introduction 

There  is  interest  in  research  on  nitramine  combustion  because  among  energetic  materials, 
nitramines  possess  many  desirable  qualities  such  as  high  specific  impulse  and  stable  combus¬ 
tion  in  propulsive  systems  where  they  are  used  and  low  levels  of  pollutants  and  smoke  formed 
during  combustion.  These  attributes  have  prompted  many  experimental  and  theoretical  stud¬ 
ies  on  the  combustion  of  nitramines.  The  physical  and  chemical  processes  that  govern  the 
ignition  and  combustion  of  nitramine  propellants  are  complex  and  involve  decomposition  of 
the  propellant  into  several  compounds  and  combustion  of  these  compounds  in  stages.  Kubota 
(1981)  observed  that  the  nitramine  particles  melt  and  diffuse  into  a  melting  binder  on  the 
burning  surface  of  the  propellant  and  generate  a  homogeneous  mixture  of  reactive  species. 
This  homogeneous  mixture  then  oxidizes  in  the  gas  phase.  Melius  (1988)  identified,  in  the  gas 
phase,  a  primary  flame  zone  composed  of  a  very  endothermic  pyrolysis  region  followed  by  a 
dark  zone  and  a  luminous  secondary  reaction  zone.  Previous  numerical  studies  using  detailed 
gas-phase  chemistry  have  demonstrated  the  importance  of  the  secondary  reaction  zone  in  con¬ 
trolling  the  burning  rate  of  neat  RDX  (Melius  1986,  Melius  1987,  Melius  1988,  Melius  1989). 
The  success  of  these  numerical  studies  illustrates  the  utility  of  fundamental  approaches  to 
clarify  the  combustion  of  nitramines. 

The  gas-phase  combustion  of  nitramines  takes  place  in  a  number  of  subsystems,  and  the 
reactions  between  CH4-NO2  and  CH2O-NO2  are  such  subsystems.  An  improved  understand¬ 
ing  of  the  chemistry  taking  place  in  these  systems  can  be  useful  in  modeling  the  combustion 
of  nitramines.  Hydrocarbon  compounds  evolve  from  decomposition  of  the  binder  and  CH4 
is  chosen  to  represent  these  hydrocarbon  compounds.  In  solid  propellants  where  the  fuel 
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and  oxidizer  are  not  premixed  on  the  molecular  level,  they  evolve  separately  and  react  in 
a  nonpremixed  configuration.  For  this  reason  numerical  calculations  are  performed  in  CH4> 
NO2  systems.  These  calculations  are  performed  using  a  detailed  chemical-kinetic  mechanism, 
a  skeletal  chemical- kinetic  mechanism,  a  10-step  reduced  chemical- kinetic  mechanism  and  a 
7-step  reduced  chemical-kinetic  mechanism.  The  skeletal  mechanism  is  deduced  from  the  de¬ 
tailed  mechanism  after  removing  reactions  which  are  found  to  have  a  negligible  influence  on 
the  flame  structure  and  critical  conditions  of  extinction.  The  reduced  mechanisms  are  derived 
from  the  skeletal  mechanism. 

As  an  initial  step  toward  understanding  gas-phase  nitramine  combustion,  results  of  numer¬ 
ical  calculations  of  the  structure  of  counterflow,  nonpremixed  methane/nitrogen  dioxide  flame 
studies  are  discussed  for  various  values  of  pressure.  The  primary  issue  addressed  in  this  work 
is  the  influence  of  strain  on  the  structure  and  extinction  of  these  flames  for  a  pressure  range 
from  1  to  10  atm.  Although  experimental  studies  on  premixed  flames  burning  CH4/NO2/O2 
mixtures  are  available  (Branch,  Sadeqi,  Alfarayedhi  &  Tiggelen  1991,  Zabarnick  1991),  similar 
experimental  studies  on  nonpremixed  CH4/NO2  flames  have  not  yet  been  conducted. 


2.2  Detailed  Chemical  Kinetic  Mechanism  and  Numerical  For¬ 
mulation 

A  detailed  chemical-kinetic  mechanism  comprising  137  elementary  reactions  among  33  species 
shown  in  Table  2.1  is  used  to  describe  the  structure  of  CH4/NO2  flames.  All  elementary  re¬ 
actions  are  presumed  to  be  reversible.  For  each  elementary  reaction  n,  Table  2.1  shows  the 
values  of  the  frequency  factor  the  temperature  exponent  6^,  and  the  activation  energy  En 
that  appear  in  the  parametric  expression  kn  =  AnT^^exp[—En/{RT]  for  the  specific  reaction 
rate  constant  kn-  Here  R  represents  the  universal  gas  constant  and  T  the  gas  temperature. 
This  mechanism  is  deduced  from  the  comprehensive  mechanism  of  Miller  &  Bowman  (1989) 
by  including  reactions  appropriate  to  CH4/NO2  systems.  The  mechanism  of  Miller  &  Bow¬ 
man  (1989)  was  developed  to  describe  the  gas  phase  reactions  of  nitrogen  compounds  that  are 
required  to  predict  the  formation  of  air  pollutants.  This  mechanism  (Miller  &  Bowman  1989) 
did  not  include  some  reactions  which  could  be  important  is  systems  where  NO2  and  NO  are 
present  in  large  quantities  such  as  those  considered  here.  To  make  up  for  this  exclusion  re¬ 
actions  132-137  taken  from  Zabarnick  (1991)  are  included  in  the  detailed  mechanism.  The 
elementary  reaction  1  taken  from  Smooke  &  Giovangigh  (1991a)  is  also  included  in  the  detailed 
mechanism  because  this  reaction  is  found  to  influence  the  critical  conditions  of  extinction  of 
methane-air  flames. 
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In  developing  the  chemical-kinetic  mechanism  reactions  which  include  species  comprising 
two  or  more  carbon  atoms  are  neglected.  This  approximation  is  consistent  with  previous  nu¬ 
merical  studies  on  the  structure  of  nonpremixed  methane-air  flames  (Smooke,  Puri  &  Seshadri 
1986,  CheUiah,  Law,  Ueda,  Smooke  Williams  1990,  Smooke  &  GiovangigU  1991a,  Smooke  & 
Giovangigli  1991b,  CheUiah,  Seshadri  &  Law  1993).  In  nonpremixed  CH4/NO2  flames  it  could 
be  assumed  that  species  comprising  two  or  more  carbon  compounds  are  formed  on  the  fuel 
side  of  the  flame  as  in  methane- air  flames,  foUowing  the  pyrolysis  of  methane  and  the  forma¬ 
tion  of  methyl  radicals.  However,  numerous  previous  studies  that  provided  valuable  insight  on 
methane-air  flames  restricted  the  chemical  kinetic  representation  to  include  only  Ci  chemistry 
(by  neglecting  reactions  with  species  comprising  two  or  more  carbon  compounds).  The  flame 
structure  and  critical  conditions  of  extinction  calculated  using  this  simplified  representation 
were  in  good  agreement  with  measurements  (Smooke  et  al.  1986,  CheUiah  et  al.  1990).  Fur¬ 
thermore  Zabarnick  (1991)  cited  the  lack  of  high-temperature  kinetic  data  for  the  reactions 
of  NO2  with  hydrocarbon  species  as  a  probable  cause  for  the  disagreement  between  modeUng 
and  experimental  results  of  CH4/NO2/O2  premixed  flames.  This  work  by  Zabarnick  (1991) 
used  the  chemical-kinetic  mechanism  of  MiUer  &  Bowman  (1989)  as  the  base  mechanism.  In¬ 
deed,  there  appears  to  be  scarce  high  temperature  data  (above  600  K)  for  the  reactions  among 
species  comprising  two  or  more  carbon  atoms  and  nitrogen  compounds. 

These  observations  suggest  that  either  the  the  chemical  reactions  taking  place  on  the  fuel 
side  of  nonpremixed  CH4/NO2  flames  is  similar  to  those  on  the  fuel  side  of  nonpremixed 
methane/air  flame,  or,  less  Ukely,  that  if  they  were  different,  the  reactions  and  their  rates  are 
not  known  with  sufficient  accuracy.  Therefore  the  chemical- kinetic  mechanism  employed  here 
includes  only  Ci  chemistry. 

Numerical  calculations  are  performed  to  determine  the  structure  and  critical  conditions  of 
extinction  of  diffusion  flames  stabiUzed  between  two  count erflowing  streams,  one  containing 
methane  and  the  other  containing  nitrogen  dioxide.  The  ambient  reactant  stream  containing 
the  fuel  is  presumed  to  flow  toward  the  stagnation  plane  from  y  =  —00,  and  the  ambient 
reactant  stream  containing  nitrogen  dioxide  is  presumed  to  flow  toward  the  stagnation  plane 
from  y  =  00.  Here  y  represents  the  spatial  coordinate  normal  to  the  surface  of  the  flame. 
Smooke  h  Giovangigli  (1991a)  have  summarized  the  conservation  equations  for  mass,  mo¬ 
mentum  and  energy,  and  the  balance  equation  for  species  which  are  used  here  to  describe 
the  flame  structure.  In  configurations  where  the  external  flow  is  inviscid  and  irrotational,  a 
parameter  reflecting  the  fluid  dynamic  environment  can  be  taken  to  be  the  strain  rate  exerted 
on  the  flame  by  the  external  oxidizer  stream  Kq  =  (—dvy/dy)^  =  jaoo  where  Vy  represents 
the  velocity  component  normal  to  the  surface  of  the  flame,  and  the  subscript  00  denotes  con¬ 
ditions  in  the  external,  ambient  N02-stream.  The  quantity  j  =  1  for  planar  flow  and  j  =  2 
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for  axisymmetric  flow  and  aoo  =  {dvxj dx)oo^  where  and  x  represent  the  velocity  component 
and  spatial  coordinate  parallel  to  the  surface  of  flame  respectively. 

The  numerical  calculations  can  be  performed  employing  either  the  “potential-flow”  or  the 
“plug-flow”  boundary  conditions.  If  the  outer  flow  is  inviscid  and  irrotational  then  potential- 
flow  boundary  conditions  must  be  used.  For  potential  flow,  on  each  side  of  the  stagnation 
plane  the  value  of  dvxjdx  is  independent  of  y,  and  at  the  ambient  nitrogen  dioxide  stream 
Vxl{a^x)  =  1.0,  and  at  the  ambient  fuel  stream  Vxl{(ioo^)  =  \/Poo/ P-oo^)  where  p  is  the  gas 
density  and  the  subscript  —  cxd  denotes  conditions  in  the  external,  ambient  fuel  stream.  Also 
for  these  potential  flow  boundary  conditions  the  stagnation  plane  is  located  at  y  =  0,  where 
Vy  =  It  is  weU  known  that  “plug-flow”  boundary  conditions  must  be  used  if  the  results 
of  numerical  calculations  are  to  be  compared  with  experimental  measurements  (Seshadri  & 
WiUiams  1978,  CheUiah  et  al.  1990,  Trees,  Brown,  Seshadri,  Smooke,  Balakrishnan,  Pitz,  Gio- 
vangigli  h  Nandula  1995,  Grudno  &  Seshadri  1995).  In  the  plug-flow  approximation  =  0 
at  the  exit  of  the  ducts  from  which  the  reactants  are  introduced  into  the  flow-fleld.  Seshadri 
Sz  WiUiams  (1978)  have  shown  that  for  plug-flow,  if  the  Reynolds  number  at  the  exit  of  the 
duct  is  sufficiently  large,  then  two  inviscid  and  rotational  flow  regions  separated  by  a  thin 
viscous  region  wiU  be  formed.  Also,  in  the  plug-flow  configuration  dv^/dx  is  a  function  of  y 
and  Seshadri  &  Williams  (1978)  have  deduced  expressions  for  calculating  the  value  of  dv^/dx 
at  the  stagnation  plane.  Since  the  focus  of  this  paper  is  to  clarify  the  influence  of  strain  on 
the  structure  and  critical  conditions  of  extinction  of  nonpremixed  CH4/NO2  flames,  and  is  not 
concerned  with  comparing  results  of  numerical  calculations  with  measurements,  calculations 
are  performed  employing  potential  flow  boundary  conditions  for  flames  stabilized  in  planar 
flows. 

Numerical  calculations  are  performed  with  the  values  of  the  mass  fractions  of  methane  and 
nitrogen  dioxide  in  the  ambient  streams  set  equal  to  unity  and  the  value  of  the  temperature 
set  equal  to  330  K.  The  value  of  aoo  is  also  specified.  Calculations  are  performed  for  values 
of  pressure  p  equal  to  1,  2,  and  10  atmospheres. 


2.3  Flame  Structure 

Figures  2.1,  2.2,  and  2.3  show  profiles  of  temperature  and  the  mole  fraction  of  major  species, 
some  intermediate  species,  and  some  radicals  calculated  using  the  detailed  chemical-kinetic 
mechanism  for  aoo  =  60  s~^  and  p  =  1  atm.  The  profiles  in  Figures  2.1,  2.2,  and  2.3  are 
plotted  as  functions  of  the  mixture  fraction  Z  =  Zc I  Zc^F-i  where  Zc  is  the  mass  fraction  of 
element  carbon  and  Zc,f  is  the  value  of  Zc  at  the  ambient  fuel  stream.  It  follows  from  this 
definition  that  Z  =  0  in  the  ambient  NO2  stream  and  Z  =  1.0  in  the  ambient  fuel  stream. 
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The  mass  fraction  of  carbon  is  calculated  from  the  relation  Zc  —  IWi  where  Wc 

and  Wi  are  respectively  the  molecular  weight  of  carbon  and  species  i,  Yi  is  the  mass  fraction 
of  species  z,  ai  is  the  number  of  carbon  atoms  in  a  molecule  of  species  and  N  is  the  total 
number  of  species. 

In  Fig.  2.1  the  profiles  of  the  temperature  T  of  the  reactants  CH4  and  NO2  and  the  ma¬ 
jor  products  CO2,  H2O  and  N2  are  plotted.  The  temperature  profile  in  Fig.  2.1  shows  that 
chemical  reactions  are  taking  place  in  a  broad  zone  in  comparison  to  those  in  nonpremixed 
methane-air  flames  (Smooke  et  al.  1986).  The  fuel  and  nitrogen  dioxide  are  consumed  in 
two  diiferent  regions  of  the  flame.  The  peak  values  of  the  products  CO2,  H2O,  and  N2  are 
attained  between  the  regions  where  CH4  and  NO2  are  consumed.  Figure  2.2  shows  profiles 
of  the  intermediate  species  O2,  H2,  CO,  NO,  and  HCN.  The  peak  values  of  the  intermediate 
species  H2,  and  CO  are  attained  at  the  region  where  CH4  is  consumed  and  the  peak  values 
of  O2  and  NO  are  attained  at  the  region  where  NO2  is  consumed.  Figure  2.2  shows  the  peak 
value  of  HCN  is  reached  at  the  position  where  CH4  is  consumed.  Figure  2.3  shows  profiles  of 
the  radicals  OH,  H,  0,  and  N.  After  comparing  the  profiles  in  Figures  2.1,  2.2,  and  2.3  it  can 
be  observed  that  the  peak  values  of  the  radicals  OH,  H,  0  and  N  are  reached  between  the 
regions  where  CH4  and  NO2  are  consumed.  Also  these  radicals  are  consumed  in  the  regions 
where  CH4  and  NO2  are  consumed. 

The  profiles  in  Figures  2.1,  2.2,  and  2.3  is  used  to  divide  the  zone  where  chemical  reactions 
are  taking  place  into  three  regions — a  fuel- consumption  region,  a  N02“Consumption  region  and 
a  product-formation  region.  In  the  fuel- consumption  region  the  predominant  reactions  taking 
place  are  those  between  the  fuel  and  the  radicals  to  form  the  intermediate  compounds  H2 
and  CO.  In  the  N02-consumption  region  the  predominant  reactions  are  between  the  NO2  and 
radicals  to  form  the  intermediate  compounds  NO  and  O2.  In  the  product-formation  region  the 
products  CO2,  H2O  and  N2  are  formed  from  the  intermediate  species.  The  radicals  are  also 
formed  in  the  product-formation  region.  The  asymptotic  structure  of  the  fuel- consumption 
region  and  the  product-formation  region  appears  to  be  similar  to  the  structure  of  these  zones 
in  methane-air  flames  (Seshadri  &  Williams  1994).  These  observations  are  used  to  derive  the 
reduced  chemical-kinetic  mechanism. 

2.4  Skeletal  and  Reduced  Chemical  Kinetic  Mechanisms 

A  skeletal  chemical-kinetic  mechanism  is  derived  from  the  detailed  chemical-kinetic  mecha¬ 
nism  after  removing  species  and  reactions  which  have  a  negligible  influence  on  the  critical 
conditions  of  extinction.  AH  elementary  reactions  in  which  the  species  CH3O,  H2O2,  NH, 
H2CN,  NNH,  C2N2,  HOCN,  and  HNCO  participate  are  removed.  In  addition  the  elementary 
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steps  132-137  are  also  removed.  The  skeletal  mechanism  comprises  reactions  1-92  shown  in 
Table  2.1.  Twenty-five  species  are  present  in  the  skeletal  mechanism.  If  reactions  in  which 
HCN  and  HCNO  participate  are  removed,  then  the  value  of  of  extinction  strain  rate  is  approx¬ 
imately  5  times  lower  than  its  calcualted  value  with  these  reactions  included  in  the  mechanism. 


Procedures  for  deriving  reduced  chemical-kinetic  mechanisms  are  described  in  detail  else¬ 
where  (Peters  1991).  Since  the  skeletal  mechanism  is  made  up  of  25  reactive  species  and  four 
chemical  elements,  twenty-one  independent  species  balance  equations  can  be  written  (Seshadri 
Sz  Williams  1994).  Therefore  the  skeletal  mechanism  is  equivalent  to  a  reduced  chemical-kinetic 
mechanism  having  21  overall  steps.  To  identify  species  for  which  steady-state  approximations 
can  be  introduced,  the  magnitudes  of  the  various  terms  in  the  balance  equation  for  each  in¬ 
termediate  species  were  plotted  as  a  function  of  Z  using  the  results  of  numerical  calculations 
with  the  skeletal  mechanism.  Steady-state  approximations  are  introduced  for  those  species 
for  which  the  absolute  values  of  the  chemical  production  rates  and  destruction  rates  are  much 
larger  than  the  absolute  values  of  the  other  terms  in  the  balance  equation  (Williams  1985). 
Using  this  procedure,  the  results  of  the  numerical  calculations  show  that  it  is  reasonable  to 
introduce  steady-state  approximation  for  the  eleven  species  HO2,  CH2O,  HCO,  CH2,  CH,  CN, 
N,  HCNO,  HNO,  N2O,  and  NCO  appearing  in  the  skeletal  chemical-kinetic  mechanism.  The 
reduced  chemical-kinetic  mechanism  now  comprises  ten  overall  steps  which  can  be  written  as 


CH4  +  H  ^ 
CH3  +  O  ^ 
H  +  H  +  M  ^ 
OH  +  CO  ^ 
H  +  O2  ^ 
H  +  NO2  ^ 
NO  +  NO  ^ 
H2  +  OH  ^ 
H2  +  0  ^ 

HCN  +  NO  +  O2  +  H  ^ 


CH3  +  H2 

I 

H  +  H2  +  CO, 

II 

H2  +  M, 

III 

H  +  CO2, 

IV 

OH  +  0, 

V 

OH  +  NO, 

VI 

O2  +  N2, 

VII 

H  +  H2O, 

VIII 

H  +  OH, 

IX 

H2  +  2  0  +  CO  +  N2, 

X 
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reactions  n  =  1-92  are 


Wj  =  —W1  +  W2  +  W3  +  W4, 

WJI  =  W5  +  ^^^49  +  W59  +  W72, 

Will  =  Wi-  Ws  +  tnio  +  W12  +  Wi3  +  Wi4  +  Wi5  +  U)16+ 

W24  +  W29  +  W30  +  W31  +  1032  +  ^«33  +  ^^^34  +  W35- 

W4I  +  W52  +  ^55  +  4056  +  Wq3  +  4067  +  ^^^68  +  ^70  + 

4071  +  W7S  +  W83  +  4084  +  ^^^85  +  ^1^90  +  ^91 

WIV  —  4014  +  ^i^l6  +  WX7  +  4018  +  ^<^19  “  W37  -  W50  +  W53  +  4O55  -  4065, 

WV  =  Wis  +  Wl9  +  W20  -  W22  +  If 26  -  'if 35  “  if 37  +  ^^38  "  'if40- 

4044  +  4045  +  4046  +  7048  +  W33  +  4O54  +  4O55  +  4O56  +  4O57+ 

7f58  +  7f62  -  'if 67  “  "if 70  “  'if 71  +  'if 72  +  "if 74  +  7f75  +  ^7^- 


WVI  = 

wvn  = 

wvin  = 


WIX  = 


Wx  = 


W77  -  4078  -  ■if  81  -  7f82  -  'if 83  -  'if 84  -  "if  85  “  If  86, 

-4038  +  7f39  +  7f40  +  7f41  +  7f82, 

-4037  -  4043  -  4044  +  4045  +  4046  +  W48- 
W71  +  W72  —  4084  ~  '^85  ~  ■if86 

4O4  +  Wq  +  40io  +  4021  +  'if25  "I"  “^28  +  +  'if50  +  'if51  +  'if52  + 

W53  +  7f54  +  W55  +  24056  +  W57  +  4058  “  ■if59  +  7f60  +  7f61  +  7062+ 

7f63  +  'if64  +  7065  +  7074  +  24075  +  '^76  ~  7O77  +  W79  +  4O90, 

702  +  409  +  7013  +  7014  +  'ifl5  +  70i6  +  7020  +  7023  +  W24  -  7025“ 
7028  +  7034  +  4035  -  4036  “  7041  “  7O45  -  4O46  -  W49  +  4O51  +  4O52+ 

7f53  +  7055  +  4056  +  4057  +  4058  “  7O59  +  4063  +  7067  +  7068  +  7070+ 

7071  —  7072  +  4076  +  7078  +  ^79  +  7083  +  7084  +  7085, 

-7045  -  7046  -  'if48  -  7072  “  7073  +  7075  +  7076  “  7077  +  7079- 


(2.1) 


To  further  simplify  the  mechanism,  steady-state  assumptions  are  introduced  for  the  species 
OH,  0,  and  HCN  which  gives  seven  overall  steps  which  can  be  written  as 


CH4  +  H 
CH3  +  H  +  H2O 
H  +  H  +  M 
H2O  +  CO 
3  H2  +  O2 
H2  +  NO2 
NO  +  NO 


CH3  +  H2 

P 

3  H2  +  CO, 

IP 

H2  +  M, 

IIP 

H2  +  CO2, 

IV* 

2  H  +  2  H2O, 

V* 

H2O  +  NO, 

VP 

O2  +  N2, 

VIP 

The  rates  of  the  overall  steps  of  the  seven  step  mechanism  Wk,  k  =  P-VIP  are  respectively 
identical  to  the  rates  for  the  overall  steps  I-VII  for  the  ten-step  mechanism  shown  in  Eq.  2.1. 
The  structure  and  critical  conditions  of  flame  extinction  are  calculated  numerically  using  the 
skeletal  mechanism  and  the  reduced  chemical-kinetic  mechanisms.  The  results  are  discussed 


in  the  next  section. 
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2.5  Results  and  Discussions 


Figures  2.4,  2.5,  2.6,  2.7,  and  2.8  show  profiles  of  temperature  and  the  mole  fraction  of  re* 
actants,  products,  intermediate  species  and  radicals  calculated  using  the  detailed  chemical- 
kinetic  mechanism  (solid  line),  skeletal  chemical-kinetic  mechanism  (dashed  line),  reduced 
10-step  mechanism  (chain  line),  and  reduced  7-step  mechanism  (dotted  line)  for  aoo  =  350  s“^ 
and  p  =  1  atm.  Figures  2.4,  and  2.5  show  that  the  profiles  of  T,  CH4,  NO2,  CO25  H2O  and  N2 
calculated  using  the  skeletal  mechanism  and  the  reduced  10-step  mechanism  agree  well  with 
those  calculated  using  the  detailed  mechanism.  The  peak  value  of  the  temperature  profile 
calculated  using  the  reduced  7-step  mechanism  is  lower  than  those  calculated  using  the  other 
mechanisms  and  it  is  reached  at  a  smaller  value  of  Z.  The  profile  of  the  mole  fraction  of  the 
reactant  NO2  calculated  using  the  7-step  mechanism  show  that  this  compound  is  consumed 
at  a  lower  value  of  Z  in  comparison  to  those  calculated  using  the  other  mechanisms.  Also  the 
profile  of  N2  calculated  using  the  reduced  7-step  mechanism  deviates  from  those  calculated 
using  the  other  mechanisms.  This  is  attributed  to  inaccuracies  associated  with  steady-state 
approximations  for  OH,  0,  and  HCN. 

Figure  2.6  shows  the  profiles  of  the  intermediate  species  NO,  O2,  CO  and  H2  calculated 
using  the  skeletal  mechanism  and  the  reduced  10-step  mechanism  to  agree  well  with  those  cal¬ 
culated  using  the  detailed  mechanism.  The  profiles  of  CO  and  H2  calculated  using  the  reduced 
7-step  mechanism  agrees  reasonably  well  with  those  calculated  using  the  other  mechanisms. 
The  peak  values  of  NO  and  O2  calculated  using  the  reduced  7-step  mechanism  is  higher  than 
those  calculated  using  the  other  mechanisms  and  it  is  reached  at  a  smaller  value  of  Z,  The 
shift  in  the  locations  of  the  peak  values  of  NO  and  O2  are  consistent  with  the  shift  in  the 
location  of  the  region  where  NO2  is  consumed  as  shown  in  Fig.  2.4.  Figures  2.7,  and  2.8  show 
profiles  of  the  radicals  H,  OH  and  0.  The  peak  values  of  these  radicals  calculated  using  the 
skeletal  and  reduced  10-step  mechanism  are  in  reasonably  good  agreement.  The  peak  values  of 
the  radicals  calculated  using  the  reduced  7-step  mechanism  are  significantly  higher  than  those 
calculated  using  the  skeletal  mechanism  and  the  reduced  10-step  mechanism.  The  qualitative 
features  of  the  flame  structure  shown  in  the  profiles  in  Figs.  2.4,  2.5,  2.6,  2.7,  and  2.8  are 
similar  to  those  shown  in  Figs.  2.1,  2.2,  and  2.3.  These  figures  show  the  chemical  reactions 
to  take  place  in  three  regions — the  fuel- consumption  region,  the  N02-consumption  region  and 
the  product-formation  region. 

The  detailed  skeletal  and  reduced  chemical-kinetic  mechanisms  are  used  to  calculate  the 
critical  conditions  of  flame  extinction.  The  flame  structure  is  calculated  for  various  values 
of  the  strain  rate.  The  critical  value  of  the  strain  rate  at  extinction  is  defined  as  that  value 
beyond  which  a  converged  solution  is  not  obtained.  Table  2.2  shows  the  values  of  the  critical 
conditions  of  extinction  calculated  using  the  various  mechanisms  for  p  =  1,  2,  and  10  atm.  The 
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values  of  a^o  calculated  using  the  various  mechanisms  are  in  good  agreement.  All  mechanisms 
show  the  critical  strain  rate  at  extinction  to  increase  with  increasing  values  of  p. 

In  deducing  the  reduced  7-step  mechanism  from  the  reduced  10-step  mechanism  steady- 
state  approximations  were  introduced  for  the  species  0,  OH,  and  HCN.  The  peak  values  of 
the  flame  temperature  and  the  values  of  the  critical  conditions  of  flame  extinction  are  found 
not  to  be  influenced  by  introducing  steady-state  approximation  for  0.  The  peak  values  of  the 
flame  temperature  and  the  values  of  the  critical  conditions  of  flame  extinction  are  found  to 
increase  if  steady-state  approximation  is  introduced  for  OH  and  the  values  of  these  quantities 
are  found  to  decrease  if  steady-state  approximation  is  introduced  for  HCN. 

2.6  Summary  and  Conclusions 

Numerical  calculations  are  performed  using  a  detailed,  skeletal,  reduced  10-step  and  a  re¬ 
duced  7-step  mechanism  to  calculate  the  structure  and  critical  conditions  of  extinction  of 
nonpremixed  CH4/NO2  flames.  The  flame  structure  calculated  using  the  skeletal  mechanism 
and  reduced  10-step  mechanism  agree  well  with  those  calculated  using  the  detailed  mechanism. 
The  peak  values  of  the  temperature  and  mole  fraction  of  the  radicals  calculated  using  the  re¬ 
duced  7-step  mechanism  are  reached  shghtly  at  a  smaller  value  of  Z  than  those  calculated 
using  the  other  mechanisms.  The  critical  conditions  of  flame  extinction  calculated  using  the 
various  mechanisms  agree  reasonably  weU  with  each  other.  The  differences  in  the  values  of  the 
the  strain  rate  at  extinction  calculated  using  the  detailed  mechanism  and  skeletal  mechanism 
are  less  than  4%,  those  between  the  skeletal  mechanism  and  the  reduced  10-step  mechanism 
and  between  the  the  reduced  10-step  and  7-step  mechanisms  are  less  than  12  %. 
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Table  2.1:  Detailed  and  Skeletal  Chemical- Kinetic  Mechanisms  and  As¬ 
sociated  Rate  Constants. 


Reactions 

An 

b„ 

En 

1. 

CH3+H+M=CH4+M 

/H2  2.0/co  2.O/CO2  3.O/H2O  5.0/ 

2.11E+14 
ko  6.26E+23 

0.0 

-1.80 

0.0 

0.0 

2. 

CH4  +  02=CH3  +  H02 

7.90E+13 

0.0 

56000.0 

3. 

CH4-fH=CH3+H2 

2.20E+04 

3,0 

8750.0 

4. 

CH4+0H=CH3+H20 

1.60E+06 

2.1 

2460.0 

5. 

CH3+0=CH20+H 

8.00E+13 

0.0 

0.0 

6, 

CH20+0H=HC0+H20 

3.43E+09 

1.2 

-447.0 

7. 

CH20+H=:HC0-i-H2 

2.19E+08 

1.8 

3000.0 

8. 

CH20-fM=HC0-f-H+M 

3.31E+16 

0.0 

81000.0 

9. 

CH20+0=HC0+0H 

1.80E+13 

0.0 

3080.0 

10. 

HC0+0H=H204-C0 

l.OOE+14 

0.0 

0.0 

11, 

HCO+M=H+CO+M 

/CO  I.9/H2  I.9/CH4  2.8/CO2  3.O/H2O  5.0/ 

2.50E+14 

0.0 

16802.0 

12. 

HC0+H=C0+H2 

1.19E+13 

0.2 

0.0 

13. 

HCO+0=CO-hOH 

3.00E+13 

0.0 

0.0 

14. 

HC0+0=C02+H 

3.00E+13 

0.0 

0.0 

15. 

HC0+02=H02+C0 

3.30E+13 

-0.4 

0.0 

16. 

C0+0+M=C02+M 

6.17E+14 

0.0 

3000.0 

17. 

C0+0H=C02+H 

1.51E+07 

1.3 

-758.0 

18. 

C0+02=C02+0 

1.60E+13 

0,0 

41000.0 

19. 

H02-hC0=C02+0H 

5.80E+13 

0.0 

22934.0 

20. 

H2“j”02=20H 

1.70E+13 

0.0 

47780.0 

21. 

0H+H2=H204-H 

1.17E+09 

1.3 

3626.0 

22. 

0+0H=02+H 

4.00E+14 

-0.5 

0.0 

23. 

0-fH2=0H+H 

5.06E+04 

2.7 

6290.0 

24. 

H-t-02+M=H02+M 

/CO  2.I/H2  2.9/CO2  4.2/H2O  I.86/N2  1.3/ 

3.61E+17 

-0.7 

0.0 

25. 

OH+H02=H20-h02 

7.50E+12 

0.0 

0.0 

26. 

H+H02=20H 

1.40E+14 

0.0 

1073.0 

27. 

0-f-H02~02”l“0Il 

1.40E+13 

0.0 

1073.0 

28. 

20H=0-hH20 

6.00E+08 

1.3 

0.0 

29. 

2H+M=H2+M 
/H2  O.O/CO2  O.O/H2O  0.0/ 

l.OOE+18 

-1.0 

0.0 

30. 

2H-fH2=2H2 

9.20E+16 

-0.6 

0.0 

31. 

2H4-H20=H2+H20 

6.00E+19 

-1.2 

0,0 

32. 

2H'-hC02=H2“hC02 

5.49E+20 

-2,0 

0.0 

33. 

H-fOH+M=H20-{-M 
/H2O  5.0/ 

1.60E+22 

-2.0 

0.0 

34. 

H+0-fM=OH-hM 
/H2O  5.0/ 

6.20E+16 

-0.6 

0.0 

35. 

20+M=02+M 

1.89E+13 

0.0 

-1788.0 

36. 

H-fH02=H2+02 

1.25E+13 

0.0 

0.0 

37. 

C02+N=N0+C0 

1.90E+11 

0.0 

3400.0 

38. 

H02+N0=:N02+0H 

2.11E+12 

0.0 

-479.0 

39. 

N02-f-H=NO+OH 

3.50E+14 

0.0 

1500.0 

40. 

N02-f0==N0-^02 

l.OOE+13 

0.0 

600.0 

41. 

N02+M=:N0  +  0  +  M 

l.lOE+16 

0.0 

66000.0 

42. 

N+N0=N2+0 

3.27E+12 

0.3 

0.0 

43, 

N+02=N0+0 

6.40E+09 

1.0 

6280.0 

44. 

N+OH=NO+H 

3.80E+13 

0.0 

0.0 

45. 

CH+N2=HCN+N 

2.00E+11 

0.0 

13600.0 

46. 

CH+NO=HCN+0 

l.lOE+14 

0.0 

0.0 
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47. 

CH2+H=CH+H2 

l.OOE+18 

-1.6 

0.0 

48. 

CH2+N0=HCN0+H 

1.39E+12 

0.0 

-1100,0 

49. 

CH3+0H=CH2-fH20 

7.50E+06 

2.0 

5000.0 

50. 

CH2+C02=CH20+C0 

l.lOE+11 

0.0 

1000.0 

51. 

CH2+0=C0+2H 

5.00E+13 

0,0 

0.0 

52. 

CH2+0=C0+H2 

3.00E+13 

0.0 

0.0 

53. 

CH2+02=:C024-2H 

1.60E+12 

0.0 

1000.0 

54. 

CH2+02=CH20+0 

5.00E+13 

0.0 

9000.0 

55. 

CH2  +  02  =  C02+H2 

6.90E+11 

0.0 

500.0 

56, 

CH2+02=C0+H20 

1.90E+10 

0.0 

-1000.0 

57. 

CH2+02=C0+0H+H 

8.60E+10 

0.0 

-500.0 

58. 

CH2+02=HC0+0H 

4.30E+10 

0.0 

-500.0 

59, 

CH3+H=CH2+H2 

9.00E+13 

0,0 

15100.0 

60. 

CH2+0H=CH+H20 

1.13E+07 

2.0 

3000.0 

61. 

CH2-f0H=CH20+H 

2.50E+13 

0.0 

0.0 

62. 

CH+02=HC0+0 

3.30E+13 

0.0 

0.0 

63. 

CH+0=CO+H 

5.70E+13 

0.0 

0.0 

64. 

CH+OH=HCO+H 

3.00E+13 

0,0 

0.0 

65. 

CH+C02=HC0+C0 

3.40E+12 

0.0 

690.0 

66. 

CH+H20=CH20+H 

4.57E+14 

-0.8 

0.0 

67. 

N20+0H=N2+H02 

2.00E+12 

0.0 

10000.0 

68. 

N20+H=N2H-0H 

7.60E+13 

0.0 

15200.0 

69. 

N20+M=N2+0+M 

1.60E+14 

0.0 

51600.0 

70. 

N20-}-0=N2H“02 

l.OOE+14 

0.0 

28200.0 

71. 

N20+0=2N0 

l.OOE+14 

0.0 

28200.0 

72. 

CH3-fNO=HCN+H20 

l.OOE+11 

0.0 

15000.0 

73. 

CH2-l-N=HCN-fH 

5.00E+13 

0.0 

0.0 

74. 

CH+N=CN-fH 

1.30E+13 

0.0 

0.0 

75. 

HCN+0H=CNH-H20 

1.45E+13 

0.0 

10929.0 

76. 

HCN+0=CN+0H 

2.70E+09 

1.6 

29200.0 

77. 

CN+H2=HCN+H 

2.95E+05 

2.5 

2237.0 

78. 

CN-fO=rCO+N 

1.80E+13 

0.0 

0.0 

79. 

HCN+0=NCO+H 

1.38E+04 

2.6 

4980.0 

80. 

CN-f02=NC0+0 

5.60E+12 

0.0 

0.0 

81. 

CN+OH=NCO+H 

6.00E+13 

0.0 

0.0 

82. 

CN+N02=NC0+N0 

3.00E+13 

0.0 

0.0 

83. 

CN+N20=NC0+N2 

l.OOE+13 

0.0 

0.0 

84. 

NC04-0=NO+CO 

2.00E+13 

0.0 

0.0 

85. 

NC0+N=N2+C0 

2.00E+13 

0.0 

0.0 

86. 

NCO+OH=NO+CO+H 

l.OOE+13 

0.0 

0.0 

87. 

NCO+M=N+CO+M 

3.10E+16 

-0,5 

48000.0 

88. 

NC0+N0=N20+C0 

l.OOE+13 

0.0 

-390.0 

89. 

HNO+M=H+NO+M 
/H2O  IO.O/O2  2.O/N2  2.O/H2  2.0/ 

1.50E+16 

0.0 

48680.0 

90. 

HN0H-0H=N0+H20 

3.60E+13 

0.0 

0.0 

91. 

HNO+HrrHs+NO 

5.00E+12 

0.0 

0.0 

92, 

HCNO+H=HCN+OH 

l.OOE+14 

0.0 

12000.0 

93. 

CH4+H02  =  CH3+H202 

1.80E+11 

0.0 

18700.0 

94. 

CH3+H02=CH30-h0H 

2.00E+13 

0.0 

0.0 

95. 

CH3+02=CH30-h0 

2.05E+19 

-1.6 

29229.0 

96. 

CH30+H=CH3+0H 

l.OOE+14 

0.0 

0.0 

97. 

CHsO+M^CHsO-hH-hM 

l.OOE+14 

0.0 

25000.0 

98. 

CH304-H=CH20-fH2 

2.00E+13 

0.0 

0.0 

99. 

CH30-h0H=CH20-}-H20 

l.OOE+13 

0.0 

0.0 

100. 

CH30-h0=CH20-f0H 

l.OOE+13 

0.0 

0.0 

101. 

CH30+02=CH20+H02 

6.30E+10 

0.0 

2600.0 
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102.  2H02=H202+02 

103.  H202+M=20H-fM 

104.  Il202“l“H=H02”hH2 

105.  Il202+0H=H20-f"H02 

106.  CH2+N2=HCN+NH 

107.  HCN+0=NH+C0 

108.  NH+02=N0+0H 

109.  NH+NOr^NsO+H 
no.  NH-f-0H=N+H20 

111.  NH+N=N2+H 

112.  NH-fH=NH-H2 

113.  NCO+H=NH+CO 

114.  H2CN-fN=N2+CH2 

115.  H2CN+M=HCN+H+M 

116.  CH3+N0=H2CN+0H 

117.  CH3+N=H2CN+H 

118.  NNH+H=N2+H2 

119.  NNH+0H=N2+H20 

120.  NNH+0=N20+H 

121.  NH+02=HN0+0 

122.  NH+OH=:HNO+H 

123.  NNH=N2+H 

124.  NNH+N0=N2+HN0 

125.  C2N2+0=NC0+CN 

126.  C2N2+0H=H0CN+CN 

127.  OH+HCN=HOCN+H 

128.  OH+HCN=HNCO+H 

129.  HOCN+H=HNCO+H 

130.  NC0+H2=HNC0+H 

131.  CN+HCN=C2N2+H 

132.  N02+N02==N0+N0+02 

133.  N02+N0=:N20+02 

134.  N0+M=N+0+M 

135.  N0+N0=N2+02 

136.  N0+HN0=N20+0H 

137.  N0+H02=HN0+02 


2.00E+12 

0.0 

0.0 

1.30E+17 

0.0 

45500.0 

1.60E+12 

0.0 

3800.0 

l.OOE+13 

0.0 

1800.0 

l.OOE+13 

0.0 

74000.0 

3.45E+03 

2.6 

4980.0 

7.60E+10 

0.0 

1530.0 

2.40E+15 

>0.8 

0.0 

5.00E+11 

0.5 

2000.0 

3.00E+13 

0.0 

0.0 

l.OOE+14 

0.0 

0.0 

5.00E+13 

0.0 

0.0 

2.00E+13 

0.0 

0.0 

3.00E+14 

0.0 

22000.0 

l.OOE+11 

0.0 

15000.0 

3.00E+13 

0.0 

0.0 

l.OOE+14 

0.0 

0.0 

5.00E+13 

0.0 

0.0 

l.OOE+14 

0.0 

0.0 

l.OOE+13 

0.0 

12000.0 

2.00E+13 

0,0 

0.0 

l.OOE+14 

0.0 

0.0 

5.00E+13 

0.0 

0.0 

4.57E+12 

0.0 

8880.0 

1.86E+11 

0.0 

2900.0 

5.85E+04 

2.4 

12500.0 

1.98E-03 

4.0 

1000.0 

l.OOE+13 

0.0 

0.0 

8.58E+12 

0.0 

9000.0 

2.00E+13 

0.0 

0,0 

2.00E+12 

0.0 

26825.0 

l.OOE+12 

0.0 

60000.0 

4.00E+20 

-1.5 

150000.0 

1.30E+14 

0.0 

75630.0 

2.00E+12 

0.0 

26000.0 

2.00E+11 

0,0 

2000.0 

Table  2.2:  Values  of  the  critical  strain  rate  at  extinction  Uo©  in  s  ^  calculated  using  the  various 
mechanisms. 


Mechanism 

Extinction  strain  rate 

latm. 

2atm. 

10  atm. 

Detailed 

522 

960 

3810 

Skeletal 

502 

935 

3700 

10-step 

545 

1040 

3955 

7-step 

580 

1170 

4260 
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Mole  Fraction 


Figure  2.1:  Profiles  of  temperature  T  and  mole-fractions  of  major  species — CH4,  NO2,  N25 
H2O,  and  CO25  calculated  using  the  detailed  chemical-kinetic  mechanism  for  =  60  s”^  and 
p  =  1  atm. 
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Temperature  [K] 


Figure  2.2:  Profiles  of  temperature  T  and  mole-fractions  of  CH4  and  NO2  and  the  intermediate 
species — O2,  H2,  CO,  NO,  and  HCN,  calculated  using  the  detailed  chemical-kinetic  mechanism 
for  Ooo  =  60  s~^  and  p  =  1  atm. 
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Temperature  [K] 


Mole  Fraction 


Figure  2.3:  Profiles  of  temperature  T  and  mole-fractions  of  the  radicals — OH,  H,  0  and  N 
calculated  using  the  detailed  chemical- kinetic  mechanism  for  Uoo  =  60  s“^  and  p  =  1  atm. 
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Temperature  [K] 


Figure  2.4:  Profiles  of  temperature  T  and  mole- fractions  of  the  reactants  CH4,  and  NO2  cal- 
ciilated  using  the  detailed  chemical-kinetic  mechanism  (solid  line),  skeletal  chemical-kinetic 
mechanism  (dashed  line),  reduced  10-step  mechanism  (chain  line),  and  reduced  7-step  mech¬ 
anism  (dotted  line)  for  a^o  —  350  s“^  and  p  =  1  atm. 
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Mole  Fraction 


Figure  2.5:  Profiles  of  mole-fractions  of  the  products  CO2,  H2O,  and  N2  calculated  using  the 
detailed  chemical-kinetic  mechanism  (solid  line),  skeletal  chemical-kinetic  mechanism  (dashed 
line),  reduced  10-step  mechanism  (chain  line),  and  reduced  7-step  mechanism  (dotted  line)  for 
ttoo  =  350  s~^  and  p  =  1  atm. 
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Figure  2.7:  Profiles  of  mole-fractions  of  the  radical  H  calculated  using  the  detailed  chemical- 
kinetic  mechanism  (solid  line),  skeletal  chemical-kinetic  mechanism  (dashed  line),  reduced 
10-step  mechanism  (chain  line),  and  reduced  7-step  mechanism  (dotted  line)  for  a^o  =  350  s“^ 
and  n  =  1  at] 
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Figure  2.8:  Profiles  of  mole-fractions  of  the  radicals  OH  and  0  calculated  using  the  detailed 
chemical-kinetic  mechanism  (solid  line),  skeletal  chemical-kinetic  mechanism  (dashed  line), 
reduced  10-step  mechanism  (chain  line),  and  reduced  7-step  mechanism  (dotted  line)  for  Uco 
d  o  =  1 


Chapter  3 


Structures  of  CH4/NO2/O2  and 
CH2O/NO2/O2  Premixed  Flames 


3.1  Introduction 

This  chapter  is  reproduced  from  Chapter  6  of  the  Ph.D  thesis  of  Yang  (1993). 

Solid  propellants  such  as  nitramines  used  in  rockets  are  complex  compounds  contain¬ 
ing  carbon,  hydrogen,  oxygen,  and  nitrogen.  The  combustion  of  solid  propellants  is  a  very 
complicated  process  that  involves  phase  changes  and  chemical  reactions,  and  in  solid  rocket 
engines,  combustion  occurs  at  very  high  pressure  (Kuo  h  Summerfield  1984,  Becker  1988). 
These  unusual  conditions  and  complexities  make  it  difficult  to  study  the  structures  of  the 
flames  experimentally  and  theoretically.  Instead  of  directly  studying  solid  propellant  com¬ 
bustion,  a  more  simple  and  fundamental  way  is  to  study  the  combustion  of  the  species  that 
are  important  in  the  combustion  of  solid  propellants  (Bui-Pham  1992,  Sadeqi  1987,  Branch 
et  al.  1991,  Zabarnick  1991).  These  simplified  combustion  systems  should  be  able  to  main¬ 
tain  the  characteristics  of  solid  propellant  flames.  Since  CH4/NO2/O2  is  a  well  studied  flame 
and  its  reaction  mechanism  is  relatively  well  established  and,  NO2  and  CH2O  are  important 
species  in  the  combustion  of  solid  propellants,  premixed  CH4/NO2/O2  and  CH2O/NO2/O2 
flames  are  chosen  here  for  investigation. 

Branch  et  al.  (1991)  and  Sadeqi  &  Branch  (1991)  studied  CH4/NO2/O2  and  CH2O/NO2/O2 
premixed  flames  both  experimentally  and  theoretically.  The  experiments  (Branch  et  al.  1991, 
Sadeqi  &  Branch  1991)  were  done  with  laminar,  premixed,  flat  flames  under  low  pressures 

55  torr).  Stable  species  were  measured  using  quartz  microprobe  and  gas  chromatography. 
Temperatures  were  measured  by  thermocouples  and  Laser- Induced  Fluorescence.  It  was  found 
that  NO2  was  just  partially  reacted  in  the  flames,  and  pure  CH4/NO2  flames  could  not  be  stabi- 


26 


lized  but  pure  CH2O/NO2  flames  could  be  stabilized.  Zabarnick  (1991)  employed  laser-induced 
fluorescence  technique  to  study  premixed  CH4/NO2/O2  flames  at  pressure  of  55  torr.  The 
chemical-kinetic  processes  of  a  burner-stabilized  flame  were  also  modeled  (Zabarnick  1991), 
but  difficulties  were  encountered  in  the  modeling  of  a  freely  propagating  flame  with  high  start¬ 
ing  concentrations  of  NO2  (Zabarnick  1991). 

Numerical  calculations  of  CH4/NO2/O2  diffusion  and  premixed  flames  were  performed 
(Bui-Pham  1992,  Ilincic  1995)  and  it  was  found  that  CH4/NO2/O2  diffusion  flame  structures 
were  quiet  different  from  those  of  hydrocarbon  flames.  For  CH4/NO2  diffusion  flames,  fuel 
and  oxidizer  were  consumed  in  two  separated  layers.  No  solution  for  pure  CH4/NO2  premixed 
flames  could  be  obtained  (Bui-Pham  1992). 

In  the  present  work,  premixed  CH4/NO2/O2  and  CH2O/NO2/O2  flame  systems  are  in¬ 
vestigated  numerically  with  full  chemical-kinetic  mechanisms  and  reduced  chemical-kinetic 
mechanisms.  Two  different  definitions  of  equivalence  ratio  were  introduced  to  study  the  in¬ 
fluences  of  NO2  on  the  burning  velocities.  The  flame  structures  under  various  conditions  are 
studied  and  the  influence  of  NO2  on  the  structure  of  the  flame  was  explored. 


3.2  Governing  Equations  and  Formulations 

The  configuration  considered  here  for  premixed  flames  is  a  planar,  steady  state,  one-dimensional, 
laminar  flow.  The  flow  direction  is  denoted  by  x  and  the  flow  Mach-number  is  small.  The 
initial  unburned  reactants  comes  from  a:  =  ~oo  and  the  products  are  located  at  x  =  00.  There 
are  no  heat  losses  from  the  flame. 

Since  we  assume  that  the  flow  Mach-number  is  small,  so  for  this  kind  flow,  the  momen¬ 
tum  equation  degenerates  to  pressure  to  be  constant  in  the  whole  flow  field.  The  governing 
equations  of  mass,  species,  and  energy  are 

Mass  :  M  —  pu  =  constant. 

Species:  +  m,-,  i=l,  2,  K, 

Energy:  Mc,f  =  ^(AS)  -  E“=i -  E”=i  ■ 

The  system  is  closed  by  the  ideal  gas  law 

p  =  pW/RT,  (3.2) 

In  the  equations,  x  denotes  the  spatial  coordinate  fixed  to  the  flame,  p  is  the  density;  u 
the  flow  velocity;  M  the  mass  flow  rate  per  unit  area.  Yi  and  Vi  are  the  mass  fraction  and  dif- 


(3.1) 
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fusion  velocity  of  species  i  respectively,  ihi  denotes  the  mass  production  rate  per  unit  volume 
of  species  i;  Cp  the  constant  pressure  heat  capacity  of  ith  species;  T  the  temperature;  A  the 
heat  conductivity  of  the  mixture;  hi  the  specific  enthalpy  of  species  z;  W  the  mean  molecular 
weight  of  the  mixture;  p  the  pressure;  and  R  is  the  universal  gas  constant. 

Boundary  conditions  are:  at  a:  — oo,  T  =  T^oo^  Yi  =  ^.nd  all  gradients  vanish 

both  at  a:  ^  —oo  and  x  ^  oo. 


The  infinite  domain  shown  here  is  difficult  to  be  used  in  practical  numerical  calculations. 
So  it  is  truncated  to  a  finite  size  of  L.  The  new  domain  is  from  x  =  0  to  x  =  L,  An  ignition 
temperature  Ti  is  introduced  at  a:  =  0.  Below  this  temperature,  the  reactions  are  frozen.  Thus 
the  boundary  conditions  used  in  practical  numerical  calculations  are  the  following 
at  a;  =  0 


Ti, 

U-e 


A  dT 
dx 


=  T^c 


x=0 


(3.3) 


Sit  X  =  L 


=  0, 


(  ___ 


)x=L 


(3.4) 


One  extra  boundary  condition  shown  here  is  used  to  get  the  flow  rate  M.  For  CH4/NO2/O2 
premixed  flames,  the  following  conditions  are  given 


T  =  r_oo  =  298  K, 

^CU^, — 009  ^T'^02  ^N02, — 00,  ^02  ^^2, — CO, 

Xi-00  =  0,  i  ^  CHi,  NO2,  O2, 
and  for  CH2O/NO2/O2  flames, 


(3.5) 


T  =  T_oo  =  423  ir, 

i/4, —  00?  ^N02  ^N02i  —  <X>^  ^02  -^02, “00? 

X,  _oo  =  0,  NO2,  O2. 


(3.6) 


In  Eqs.  3.5  and  3.6,  Xi  represents  the  molar  fraction  of  species  z.  For  CH2O/NO2/O2 
systems,  T-oo  =  423  K  instead  of  298  K  was  chosen  to  be  the  initial  temperature.  Since  in 
experiments,  CH2O  can  not  be  stored  as  gas  in  a  cylinder  because  of  its  reactivity,  particularly 
with  itself  to  form  solid  polymer,  so  a  gaseous  CH2O  generator  at  certain  temperature  (say 
423  K)  is  needed  to  supply  a  continuous  stream  of  gaseous  CH2O  (Branch  et  al.  1991). 
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In  order  to  define  the  equivalence  ratio  (/>,  we  assume  that  one-step  reactions  occur  accord¬ 
ing  to 

CH4  +  NO2  “h  O2  — ^  CO2  +  N2  +  2/^20, 

CH2O  +  \{N02  +  O2)  CO2  +  \N2  +  H2O. 

The  equivalence  ratio  is  defined  as 

.  [YfIYq] 

^  [YFlYo\t' 

here  Yi  is  the  mass  fraction  of  species  i,  the  subscripts  F  and  0  represent  fuel  and  oxidizer 
respectively,  and  the  subscript  st  represents  the  stoichiometric  condition. 

As  they  were  used  before  (Bui-Pham  1992),  two  different  systems  are  defined  here;  1)  CH4 
for  CH4/NO2/O2  flames  or  CH2O  for  CH2O/NO2/O2  flames  as  the  fuel,  and  NO2+O2  as  the 
oxidizer.  The  ratio  of  the  initial  molar  fractions  of  NO2  and  O2  is  kept  to  be  unity.  2)  CH4  + 
NO2  for  CH4/NO2/O2  flames  or,  CH2O+INO2  for  CH2O/NO2/O2  flames,  as  the  fuel,  and 
O2  as  the  oxidizer.  In  the  second  definition,  the  ratio  of  the  initial  molar  fractions  of  CH4  and 
NO2  is  kept  as  unity  for  CH4/NO2/O2  flames,  and  the  ratio  of  the  initial  molar  fractions  of 
CII2O  and  NO2  is  kept  as  two  for  CH2O/NO2/O2  flames. 

The  chemical-kinetic  mechanism  of  CH4/NO2/O2  flame  systems  is  shown  in  Table  3.1.  This 
mechanism  was  used  before  in  (Bui-Pham  1992)  for  both  premixed  and  diffusion  CH4/NO2/O2 
flames.  Table  3.2  shows  the  chemical-kinetic  mechanism  of  CH2O/NO2/O2  flame  systems. 

3.3  Results  and  Discussions 
3.3.1  CH4/NO2/O2  Flames 

The  calculations  were  done  for  the  initial  temperatures  as  specified  in  the  previous  section 
and  at  pressure  being  one  atmosphere.  Figure  3.1  shows  the  flame  temperature  as  a  function 
of  equivalence  ratio.  For  the  system  corresponding  to  the  first  definition  of  the  flame 
temperature  first  increases  with  the  increasing  value  of  (j)  and  attains  a  maximum  at  ^  =  1.1, 
then  it  decreases  quickly  with  the  increasing  value  of  (j).  For  the  system  corresponding  to  the 
second  definition  of  (j)^  after  the  flame  temperature  attains  a  maximum  value,  it  then  decreases 
gradually  and  finally  approaches  a  constant  value.  This  is  because  that  as  (j)  increases,  the 
oxidizer,  O2,  decreases  and  the  fuel  CH4  +  NO2  goes  to  a  constant  composition.  As  (j>  00, 

O2  vanishes,  and  the  fuel  goes  to  CH4  +  NO2  which  itself  contains  an  oxidizer  NO2. 


(3.8) 
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Table  3.1:  Chemical- Kinetic  Mechanism  of  CH4/NO2/O2  Flames 


NO. 

Reaction® 

Ai 

nt 

Ei 

1 

CF4  +  H  =  CH3  +  H2 

2.20E04 

3.0 

8750.0 

2 

CHi  +  OH^  CH3  +  n20 

1.60E06 

2.1 

2460.0 

3 

CH4  +  0  =  CHz  +  OH 

1.60E06 

2.36 

7400.0 

4 

GHz  ^■n  +  M  =  CHi  -f  M 

1.90E36 

-7.0 

9050.0 

5 

GHz  +  0  =  GH2O  +  H 

6.80E13 

0.0 

0.0 

6 

GHz  +  02  =  GHzO  +  0 

7.00E13 

0.0 

25652.0 

7 

GHz  +  02  =  GH2O  +  OH 

7.00E13 

0.0 

25652.0 

8 

OH2O  +  H  =  HOO  +  H2 

2.51E13 

0.0 

3991.0 

9 

OH2O  +  OH  =  HOO  +  H2O 

3.00E13 

0.0 

1195.0 

10 

OHzO  +  M  =  OH2O  +  H  +  M 

2.40E13 

0.0 

28812.0 

11 

OHzO  +  H  =  GHz  +  OH 

1.00E14 

0.0 

0.0 

12 

CH^^O  “h  =  C H2O  H2 

2.00E13 

0.0 

0.0 

13 

OHzO  +  02  =  OH2O  +  HO2 

6.30E10 

0.0 

2600.0 

14 

HOO  +  M  =  H  +  00  +  M 

1.60E14 

0.0 

14700.0 

15 

HOO  +  H  =  00 +  H2 

4.00E13 

0.0 

0.0 

16 

OO  +  OH  =  OO2  +  H 

1.51E07 

1.3 

-758.0 

17 

H  +  O2  =  OH  +  0 

2.00E14 

0.0 

16800.0 

18 

H2  +  0  =  0H  +  H 

1.80E10 

1.0 

8826.0 

19 

0H  +  H2  =  H2O  +  H 

1.17E09 

1.3 

3626.0 

20 

OH  +  OH  =  H2O  +  0 

6.00E08 

1.3 

0.0 

21 

H  +  O2  +  =  HO2  +  M 

2.30E18 

-0.8 

0.0 

22 

H  +  OH  +  M'>  =  H2O  +  M 

2.20E22 

-2.0 

0.0 

23 

H  +  H  +  =  H2  +  M 

1.80E18 

-1.0 

0.0 

24 

OH  +  HO2  =  H2O  +  O2 

2.00E13 

0.0 

1000.0 

25 

H  +  H02  =  ‘20H 

1.50E14 

0.0 

1004.0 

26 

H  +  HO2  —  Hz  +  O2 

2.50E13 

0.0 

700.0 

27 

HO2  +  HO2  —  H2O2  +  O2 

2.00E12 

0.0 

0.0 

28 

H2O2  +  M  =  OH  +  OH  +  M 

1.30E17 

0.0 

45500.0 

29 

H2O2  +  H  =  HO2  +  H2 

1.60E12 

0.0 

3800.0 

30 

H2O2  +  OH  =  H2O  +  HO2 

1.00E13 

0.0 

1800.0 

31 

NO2  +  N  =  2NO 

4.00E12 

0.0 

0.0 

32 

NO2  +  N  =  N2+02 

1.00E12 

0.0 

0.0 

33 

N  +  H02  =  NO  +  OH 

1.00E13 

0.0 

2000.0 

34 

NO2  +  M  =  NO  +  0  +  M 

1.10E16 

0.0 

65570.0 

(to  be  continued  on  next  page) 
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Table  3.1  (continued) 


NO. 

Reaction 

Ai 

n,- 

Ei 

35 

2NO2  =  2NO  +  O2 

2.00E12 

0.0 

26620.0 

36 

NO  +  M  =  N  +  0  +  M 

4.00E20 

-1.5 

150000.0 

37 

2NO  =  N2  +  02 

1.30E14 

0.0 

75630.0 

38 

NO2  +  CH3  =  CH3O  +  NO 

1.50E13 

0.0 

0.0 

39 

NO2  +  H  =  NO  +  OH 

3.50E14 

0.0 

1500.0 

40 

N  +  N0  =  N2  +  0 

3.27E12 

0.3 

0.0 

41 

N  +  02  =  NO  +  0 

6.40E09 

1.0 

6280.0 

42 

NO  +  H  =  N  +  OH 

2.63E14 

0.0 

50410.0 

43 

NO2  +  0  =  NO +  O2 

1.00E13 

0.0 

600.0 

44 

NO2  +  HCO  ^C02  +  H  +  N0 

1.50E12 

0.0 

-430.0 

45 

NO2  +  00  =  CO2  +  NO 

1.26E14 

0.0 

27600.0 

^  Rate  constants  are  in  the  form  ki  =  ex^[— Ei / (RT)],  Units  are  moles,  cubic  centime¬ 

ters,  seconds,  Kelvins  and  cal/mole. 

^  Third-body  efficiencies  are  H2O:  6.5,  CO2:  1.5,  CO:  0.75,  O2:  0.4,  N2:  0.4,  all  other  species 
1.0. 
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Table  3.2:  Chemical-Kinetic  Mechanism  of  CH2O/NO2/O2  Flames 


NO. 

Reaction^ 

Ai 

Ei 

1 

CH2O  -f  OH 

=  ECO  +  H2O 

3.00E13 

0.0 

1195.0 

2 

CE2O  +  H 

=  HCO  +  H2 

2.51E13 

0.0 

3991.0 

3 

CH2O  +  O2  - 

=  HCO  +  HO2 

2.00E13 

0.0 

38950.0 

4 

CH2O  4-  HO2 

=  ECO  -b  H2O2 

2.00E12 

0.0 

11660.0 

5 

CH2O  +  M  = 

HCO  +  H  +  M 

3.31E16 

0.0 

81000.0 

6 

C  H2O  +  0  = 

=HCO+OH 

1.80E13 

0.0 

3080.0 

7 

HCO  -t-  M  = 

H  +  CO +  M 

1.60E14 

0.0 

14700.0 

8 

ECO 

^co  +  n2 

4.00E13 

0.0 

0.0 

9 

CO  +  OH 

=  C02  +  H 

1.51E07 

1.3 

-758.0 

10 

H  +  O2-- 

=  0E  +  0 

2.00E14 

0.0 

16800.0 

11 

H2  +  O-- 

=  OH  +  H 

1.80E10 

1.0 

8826.0 

12 

OH  +  H2 

=  H20  +  H 

1.17E09 

1.3 

3626.0 

13 

OH  +  OH 

=  H20  +  0 

6.00E08 

1.3 

0.0 

14 

n  +  H  +  M'>  =  E2  +  M 

1.80E18 

-1.0 

0.0 

15 

H  +  OH  +  M'’  =  E2O  +  M 

2.20E22 

-2.0 

0.0 

16 

H  +  02  +  M>>  =  EO2  +  M 

2.30E18 

-0.8 

0.0 

17 

E  +  HO2  =  20  H 

1.50E14 

0.0 

1004.0 

18 

H  +  HO2 

=  H2  +  O2 

2.50E13 

0.0 

700.0 

19 

OH  “h  HO2 

=  H2O  +  O2 

2.00E13 

0.0 

1000.0 

20 

HO2  +  HO2 

=  H2O2  +  O2 

2.00E12 

0.0 

0.0 

21 

H2O2  “1“  M  = 

OE^OE-^M 

1.30E17 

0.0 

45500.0 

22 

E2O2  +  H 

=  HO2  +  H2 

1.60E12 

0.0 

3800.0 

23 

H2O2  +  OH 

=  H2O  +  HO2 

1.00E13 

0.0 

1800.0 

24 

NO2  +  H 

=NO+OH 

3.47E14 

0.0 

1470.0 

25 

NO  ^H 

=  N  +  OH 

2.63E14 

0.0 

50410.0 

26 

N  02^-0 

=  NO +  O2 

1.00E13 

0.0 

600.0 

27 

NO2  +  N  =  2NO 

4.00E12 

0.0 

0.0 

28 

NO2  +  N 

Z=1  N2  +  O2 

1.00E12 

0.0 

0.0 

29 

NO2  +  M  = 

NO+O+M 

1.10E16 

0.0 

65570.0 

30 

2N02^ 

2NO  +  O2 

2.00E12 

0.0 

26620.0 

31 

NO2  +  ECO  = 

■  CO2  +  H  +  NO 

1.50E12 

0.0 

-430.0 

32 

NO2  +  CO 

=  CO2  +  NO 

1.26E14 

0.0 

27600.0 

33 

COa  +  iV: 

=  NO  +  CO 

1.90E11 

0.0 

3400.0 

34 

N  +  no2-- 

=NO+OH 

1.00E13 

0.0 

2000.0 

(to  be  continued  on  next  page) 
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Table  3.2  (continued) 


NO. 

Reaction 

Ai 

Hi 

Ei 

35 

NO  +  M  =  N  +  0  +  M 

4.00E20 

-1.5 

150000.0 

36 

2NO  =  N2  +  02 

1.30E14 

0.0 

75630.0 

37 

N  +  NO  =  N2  +  O 

3.27E12 

0.3 

0.0 

38 

N  +  02  =  N0  +  0 

6.40E09 

1.0 

6280.0 

39 

CH2O  +  NO2  =  HCO  +  HONO 

2.94E11 

0.0 

12850.0 

40 

HCO  +  NO2  =  CO  +  HONO 

1.24E23 

-3.29 

2355.0 

41 

HCO  +  iVO  =  CO  +  HNO 

7.20E13 

-0.40 

0.0 

42 

CH2O  +  0  =  CO2  +  2ir 

3.50E05 

2.40 

1360.0 

43 

HCO  +  OH  =  H2O  +  CO 

3.00E13 

0.00 

0.0 

44 

OH  +  NO  +  M'>  =  HONO  +  M 

1.60E14 

-0.50 

0.0 

45 

HNO  +  OH  ^  H2O  +  NO 

1.30E12 

0.50 

1990.0 

46 

CH2O  +  M  =  fl'a  +  CO 

2.30E15 

0.50 

64800.0 

47 

H-\-NO  +  M^  ^  HNO  +  M 

5.40E15 

0.00 

-600.0 

48 

HONO  +  0.ff  =  H2O  +  NO2 

1.00E13 

0.00 

3300.0 

49 

HNO  -^H  =  H2-^N0 

1.30E13 

0.00 

3970.0 

50 

HNO  +  0  =  OH  -VNO 

5.00E11 

0.50 

1990.0 

51 

HCO  +  0  =  +  CO2 

3.00E13 

0.00 

0.0 

52 

HCO  +  0  =  OH  +  CO 

3.00E13 

0.00 

0.0 

53 

0  +  0  +  M^  =  02  +  M 

1.89E13 

0.00 

-1788.0 

54 

CO  +  0  +  M*-  =  CO2  +  M 

6.17E14 

0.00 

3000.0 

55 

CO  "i"  O2  —  CO2  "I"  0 

1.60E13 

0.00 

41000.0 

56 

HNO  +  0  =  NO2  +  H 

5.00E10 

0.50 

2000.0 

57 

HCO  +  HNO  =  CH2O  +  NO 

2.00E12 

0.00 

5000.0 

58 

H  +  HONO  =  H2  +  NO2 

1.00E12 

0.00 

1000.0 

®  Rate  constants  are  in  the  form  ki  =  AiT'^*ex-p[—Ei/{RT)].  Units  are  moles,  cubic  centime¬ 
ters,  seconds,  Kelvins  and  cal/mole. 

^  Third-body  efficiencies  are  H2O:  6.5,  CO2:  1.5,  CO:  0.75,  O2:  0.4,  N2:  0.4,  all  other  species 

1.0. 

^  Third-body  efficiencies  are  H2O  :  6.5,  CO2  :  1*5,  CO  :  0.75,  O2  '  0.4,  N2  :  0.4,  aU 
other  species  1.0. 

Figure  3.2  gives  the  plots  of  burning  velocities,  as  functions  of  equivalence  ratio  0. 
For  the  system  corresponding  to  the  first  definition  of  <^,  as  (j)  increases,  the  burning  velocity 
increases,  then  attains  a  maximum  value  around  (f>  =  0.9.  As  (j)  increases  further,  the  burning 
velocity  drops  quickly.  This  behavior  is  similar  to  that  of  a  usual  two  reactants  (fuel  and 
oxidizer)  premixed  flame  system.  While  for  the  system  corresponding  to  the  second  definition 
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Figure  3.1:  Variation  of  flame  temperature  as  a  function  of  equivalence  ratio  for  CH4/NO2/O2 
flame  system 


Figure  3.2:  Variation  of  burning  velocity  as  a  function  of  equivalence  ratio  for  CH4/NO2/O2 
flame  system 


34 


of  0,  there  are  some  interesting  different  phenomena.  As  (j>  increases,  the  burning  velocity 
increases,  then  attains  a  maximum  value  around  cj)  =  0.5.  As  (j)  increases  further,  the  burning 
velocity  decreases  gradually  to  a  constant  value  which  is  3.3  cm/s  corresponding  to  cf)  —  co. 
These  facts  suggest  that  NO2  behaviors  like  an  inhibitor  and,  a  pure  CH4+NO2  premixed 
flame  is  a  very  weak  flame.  The  reasons  will  be  explored  later  in  the  discussions  of  the  flame 
structures. 

The  structures  for  CH4/NO2/O2  flames  with  the  second  definition  of  equivalence  ratio 
are  shown  in  Figures  3.3  -  3.14.  The  corresponding  equivalence  ratios  are  0.3,  0.5,  0.9  and 
00.  <j>  =  0.3  represents  the  region  in  which  the  burning  velocity  increases  with  the  increas¬ 
ing  values  of  (j>  =  0.5  corresponds  to  the  maximum  burning  velocity,  ^  =  0.9  represents 
the  region  where  burning  velocity  decreases  with  the  increasing  values  of  (j),  kt  (j)  —  0.9  the 
burning  velocity  is  roughly  the  same  as  that  dit  (f)  =  0.3.  (j)  =  00  means  that  there  is  no 
O2  presenting  in  the  system,  the  reactants  are  CH4  and  NO2  with  an  initial  composition, 
—  •  -X^02,  — 00  =  1  •  1- 

For  (j)  =  0.3  and  0.5,  as  it  shows  in  Figures  3.3  and  3.6,  the  concentration  of  oxygen 
decreases  to  a  constant  value.  For  cj)  =0.9,  as  seen  in  Figure  3.9,  the  concentration  of  oxygen 
decreases  first  to  a  minimum  value  and  then  increases  gradually.  The  decreasing  trend  of 
oxygen  concentration  can  be  attributed  to  the  reactions  of  O2  with  the  radicals.  This  trend 
win  soon  be  balanced  or  even  suppressed  (in  the  case  (j)  =  0.9)  by  the  oxygen  production 
reaction 


2NO  =  N2  +  O2. 

It  can  be  seen  from  Figures  3.3,  3.6  and  3.9  that  NO2  and  CH4  are  consumed  almost  in  the 
same  region,  and  NO2  disappears  earlier  than  CH4.  The  trend  of  NO2  profile  shown  here  is  in 
consistent  with  those  of  experimental  measurements  and  numerical  simulations  in  (Zabarnick 
1991),  but  in  disagreement  with  the  experimental  results  in  (Sadeqi  1987,  Zabarnick  1991) 
in  which  only  partial  consumption  of  NO2  in  the  flames  was  observed.  This,  as  discussed  in 
(Zabarnick  1991),  is  likely  due  to  the  use  of  a  quartz  microprobe  and  gas  chromatography 
analysis.  It  is  possible  that  NO2  could  be  formed  in  the  probing  and  analysis  processes  due 
to  the  reaction 

2NO  +  02^  2AO2. 

The  almost  simultaneous  disappearances  of  NO2  and  CH4  make  premixed  CH4/NO2/O2  flames 
very  different  from  CH4/NO2/O2  diffusion  flames.  For  the  diffusion  flames,  as  it  was  observed 
in  (Bui-Pham  1992,  Hincic  1995),  the  consumptions  of  CH4  and  NO2  occur  in  two  different 
regions. 
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Figure  3.3:  Stable  species  profiles  at  ^  =  0.3  for  CH4/NO2/O2  flame  system  with  the  second 
definition  of  equivalence  ratio 
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Figure  3.4:  Radical  species  profiles  at  ^  =  0.3  for  CH4/NO2/O2  flame  system  with  the  second 
definition  of  equivalence  ratio 
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Figure  3.5:  Intermediate  species  profiles  at  (^  =  0.3  for  CH4/NO2/O2  flame  system  with  the 
second  definition  of  equivalence  ratio 


Figure  3.6:  Stable  species  profiles  at  </>  =  0.5  for  CH4/NO2/O2  flame  system  with  the  second 
definition  of  equivalence  ratio 
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Figure  3.7:  Radical  species  profiles  at  <?!>  =  0.5  for  CH4/NO2/O2  flame  system  with  the  second 
definition  of  equivalence  ratio 


Figure  3.8:  Intermediate  species  profiles  at  0  =  0.5  for  CH4/NO2/O2  flame  system  with  the 
second  definition  of  equivalence  ratio 
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Figure  3.9:  Stable  species  profiles  dX  (f>  =  0.9  for  CH4/NO2/O2  flame  system  with,  the  second 
definition  of  equivalence  ratio 


Figure  3.10:  Radical  species  profiles  at  ^  =  0.9  for  CH4/NO2/O2  flame  system  with  the  second 
definition  of  equivalence  ratio 
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Figure  3.11:  Intermediate  species  profiles  at  (^  =  0.9  for  CH4/NO2/O2  flame  system  with,  the 
second  definition  of  equivalence  ratio 


Figure  3.12:  Stable  species  profiles  at  ^  00  for  CH4/NO2/O2  flame  system  with  the  second 
definition  of  equivalence  ratio 
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Figure  3.13:  Radical  species  profiles  at  ^  >  oo  for  CH4/NO2/O2  flame  system  with  the  second 
definition  of  equivalence  ratio 
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Figure  3.14:  Intermediate  species  profiles  at  ^  >  00  for  CH4/NO2/O2  flame  system  with  the 
second  definition  of  equivalence  ratio 
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Maxim-um  values  of  H2O,  NO,  H2,  and  CO  occur  near  the  region  where  CH4  and  NO2 
disappear.  They  are  produced  from  the  consumptions  of  CH4,  O2  and  the  decomposition  of 
N02»  This  is  similar  to  methane-air  premixed  flames  where  a  fuel  consumption  layer  exists 
to  form  H2,  CO  and  H2O.  After  the  formation  of  CO,  NO,  and  H2,  NO  is  further  attacked 
by  the  radicals  to  form  N2,  CO  is  oxidized  to  form  CO2,  and  H2  is  oxidized  to  form  H2O. 
This  layer  is  similar  to  the  oxidation  layer  of  H2  and  CO  of  methane-air  flames  (Seshadri  & 
Peters  1988b,  Seshadri  &  Gottgens  1991,  Seshadri  &  Peters  1988a,  Trevino  h  Williams  1988). 

The  profiles  of  OH,  0,  and  H  in  Figures  3.4,  3.7  and  3.10  show  that  these  radicals  always 
disappear  at  the  same  point  where  NO2  disappears.  This  suggests  that  NO2  is  a  radical  de¬ 
stroyer  and,  indeed  it  is  a  stronger  radical  destroyer  than  CH4.  In  the  combustion  process, 
NO2  competes  with  CH4  for  radicals.  Since  CH4/NO2  premixed  flames  are  very  weak  flames, 
this  competition  makes  NO2  act  like  a  fuel  and  decreases  the  levels  of  radical  concentrations, 
and  makes  the  flame  dilRcult  to  propagate.  Hence,  the  presence  of  oxygen  in  the  flame  is 
necessary  to  maintain  a  certain  level  of  radical  concentration  and  keep  the  flame  propagating. 
This  conclusion  is  supported  by  the  calculations  of  pure  CH4  +  NO2  combustion  which  are 
shown  in  Figures  3.12  -  3.14.  The  maximum  values  of  the  molar  fractions  of  the  radicals  H, 
OH,  and  0  for  the  the  pure  CH4  +  NO2  flame  are  of  orders  of  10““^,  10“^  and  10“^,  while 
those  for  CH4/NO2/O2  flames  with  (f>  =  0.3,  0.5,  and  0.9  are  of  orders  of  10“^.  The  absence 
of  oxygen  decreases  the  levels  of  radical  concentrations  and  thus  makes  the  flame  harder  to 
propagate  than  the  flame  with  initial  oxygen  does.  The  flame  velocity  for  the  pure  CH4  +  NO2 
flame  is  3.3  cm/s,  while  for  the  flames  with  <j)  =  0.3,  0.5  and  0.9,  the  flame  velocities  are  116.1, 
142.1  and  117.9  cm/s  respectively.  They  are  much  larger  than  3.3  cm/s.  A  pure  CH4  +  NO2 
flame  is  a  very  weak  flame.  This  explains  the  difficulties  encountered  by  Bui-Pham  (1992)  in 
getting  a  pure  methane,  pure  nitrogen  dioxide  flame,  and  the  dilRculties  met  by  Zabarnick 
(1991)  in  calculating  a  freely  propagating  flame  with  high  starting  concentrations  of  NO2  and, 
why  Branch  et  al.  (1991)  could  not  stabilize  flames  in  rich  mixtures  of  CH4/NO2/O2  and 
stable  flames  in  lean  mixture  could  only  be  obtained  with  difficulty. 

The  profiles  of  intermediate  species  CH3,  CH3O,  CH2O,  HCO,  HO2,  and  H2O2  are  plotted 
in  Figures  3.5,  3.8,  3.11  and  3.13. 

The  structures  of  the  CH4/NO2/O2  flames  with  the  first  definition  of  equivalence  ratio  (f) 
are  essentially  the  same  as  those  discussed  above. 


42 


4) 


Figure  3.15:  Variation  of  flame  temperature  as  a  function  of  equivalence  ratio  for 
CH2O/NO2/O2  flame  system 


4) 


Figure  3.16:  Variation  of  burning  velocity  as  a  function  of  equivalence  ratio  for  CH2O/NO2/O2 
flame  system 
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Figure  3.17:  Stable  species  profiles  at  <?!>  =  1  for  CH2O/NO2/O2  flame  system  with  the  second 
definition  of  equivalence  ratio 


Figure  3.18:  Radical  species  profiles  at  ^  =  1  for  CH2O/NO2/O2  flame  system  with  the  second 
definition  of  equivalence  ratio 
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Figure  3.19:  Intermediate  species  profiles  at  0  =  1  for  CH2O/NO2/O2  flame  system  with  the 
second  definition  of  equivalence  ratio 


3.3.2  CH2O/NO2/O2  Flames 

Figures  3.15  and  3.16  show  the  temperature  and  burning  velocity  as  functions  of  equivalence 
ratio  for  CH2O/NO2/O2  systems.  They  are  very  similar  to  the  temperature  and  burning  ve- 
locity  profiles  shown  in  Figures  3.1  and  3.2  for  CH4/NO2/O2  systems.  The  maximum  value  of 
burning  velocity  of  CH2O/NO2/O2  systems  is  446  cm/s  at  (j>=  0.6  with  the  second  definition 
of  equivalence  ratio,  and  is  395.4  cm/s  at  (f>  ==1.0  with  the  first  definition  of  equivalence  ratio. 
For  the  pure  CH2O  +  NO2  flame  with  an  initial  composition,  XCH2O  •  ^7V02  =2:1,  which 
corresponds  to  <j!>  =  00,  the  burning  velocity  is  35.7  cm/s.  Although  this  is  small  compared 
with  the  maximum  value  of  446  cm/s,  but  it  is  close  to  the  burning  velocity  of  a  methane-air 
premixed  flame  near  stoichiometry.  So  we  can  stiU  expect  a  stable  flame  for  pure  CH2O  + 
NO2  systems,  as  it  was  observed  in  (Branch  et  al.  1991). 

Only  the  structures  corresponding  to  the  flame  with  equivalence  (^  =  1.0  were  plotted. 
The  profiles  are  shown  in  Figures  3.17  -  3.19.  They  actually  are  similar  to  the  structures  of 
CH4/NO2/O2  flames.  Indeed  the  similarities  in  flame  structures  result  in  the  similarities  in 
the  behaviors  of  burning  velocity  with  equivalence  ratio. 
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3.4  The  Reduced  Mechanism  for  Premixed 
CH4/NO2/O2  Flames 

In  the  previous  sections,  the  structures  of  premixed  CH4/NO2/O2  and 
CH2O/NO2/O2  flames  were  studied  using  full  chemical-kinetic  mechanisms.  In  this  section 
and  the  following  section,  the  structures  of  these  flames  wiU  be  studied  by  using  reduced 
mechanisms.  As  it  was  done  before,  a  reduced  mechanism  is  achieved  from  a  fuU  mechanism 
by  introducing  steady-state  assumptions  for  some  species.  Here  for  CH4/NO2/O2  premixed 
flames,  the  full  mechanism  is  the  one  shown  in  Table  3.1.  Steady-state  assumptions  are  made 
for  H2O2,  HCO,  HO2,  CH3O,  CH2O,  CH3,  0,  and  N  species.  In  numerical  calculations,  a 
steady-state  assumption  for  a  species  is  effected  by  putting  the  diffusion  and  convection  terms 
of  this  species  to  zero.  Since  no  convergence  was  achieved  for  flames  with  steady-state  as¬ 
sumption  for  species  H2O25  so  truncation  for  H2O2  in  reaction  27  of  Table  3.1,  which  is  2HO2 
=  H2O2  +  O2,  was  made  in  the  calculations.  Therefore,  a  7-step  reduced  mechanism  for 
CH4/NO2/O2  premixed  flames  can  be  represented  by 

I :  CH^  +  20 H  ^  2H2  +  H2O  +  CO 

II:  H2  +  OH  ^  H  +  H2O 

III:  0H-\-C0^H  +  C02 

IV:  H  +  H  ^  H2 

V:  E2  +  02^  20H 

VI:  H  +  N02^0H  +  N0 

VII :  2N0  —  O2  +  N2 


The  reaction  rates  of  the  global  reactions  I  -  VII  can  be  expressed  in  terms  of  elementary 
reaction  rates  in  Table  3.1  as 


wj  =  u;i  +  ^^2  +  ^3  '^4 

wii  =  W2  +  W3-  We  -  W7  +  wg-  wio  -  W12  -  -  ^17  + 

+  W22  +  'W24  +  ^^30  -  ^^34  “  ^^^36  -  ^40  '“^41  +  ^^^43 

Win  =  ^16  +  ^^^44  + 

WJV  =  W4  —  U^io  +  Wi3  +  W21  +  W22  +  ^23  “  ^^^27  +  ^29  +  ^30  “  ^^34 
-^^;35  -  Wse 

WV  =  'WeVVOiV  Wi7  +  W23  +  ^27  -  tt?29  ^^30  +  ^33  ^35  +  ^^^40 

+  W4I  -  W43 


(3.9) 
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WVI 


—  ^'U)33  +  W34  +  2ii;35  +  W3Q  +  lt;38  +  w^g  —  W40  —  W41  +  'W42 

+'U^43  +  ^44  +  ^^^45 
Wvil  =  -W31-  W33  +  W3Q  +  W37  -  W41  +  W42 

3.4.1  Comparisons  between  the  Full  Mechanism  and  the  Reduced  Mecha¬ 
nism 

Calculations  were  done  for  flame  systems  with  the  first  definition  of  equivalence  ratio.  The 
burning  velocity  as  a  function  of  equivalence  ratio,  at  pressure,  p  =  1  atm.  is  plotted  in  Fig. 
3.20.  The  solid  line  represents  the  results  from  the  full  mechanism  and  the  dotted  line  rep¬ 
resents  the  results  from  the  reduced  mechanism.  It  can  be  seen  that  they  agree  with  each  other. 

Structures  of  the  premixed  CH4/NO2/O2  flames  with  equivalence  ratio,  =  1,  and  pres¬ 
sure,  p  =  latm.  obtained  from  the  reduced  mechanism  are  shown  as  dotted  lines  in  Figs.  3.21 
-  3.25,  the  solid  lines  in  these  figures  represent  the  structures  of  the  same  flame  predicted  by 
the  full  mechanism.  The  temperature  profile  obtained  using  the  reduced  mechanism  agrees 
with  that  using  the  fuU  mechanism.  For  stable  species  NO25  CH4,  O2,  H2O,  NO,  CO,  CO2, 
N2,  H2,  and  the  radicals  OH,  H,  0,  the  results  predicted  by  the  reduced  mechanism  agree 
well  with  those  from  the  full  mechanism.  Species  profiles  for  CH3O,  CH2O,  and  CH3  are 
shown  in  Fig.  3.24.  The  peak  values  of  these  species  from  the  reduced  mechanism  are  roughly 
three  times  larger  than  those  from  the  fuU  mechanism.  Species  profiles  for  N,  HCO,  HO2, 
and  H2O2  are  shown  in  Fig.  3.25.  For  nitrogen  atom  N,  the  mole  fraction  calculated  from 
the  fuU  mechanism  are  always  higher  than  that  obtained  from  the  reduced  mechanism.  The 
peak  mole  fractions  of  species  HCO  and  HO2  predicted  by  the  reduced  mechanism  are  higher 
than  those  by  the  full  mechanism.  For  H2O2,  the  peak  mole  fraction  obtained  from  the  fuU 
mechanism  is  of  two  orders  larger  than  that  obtained  from  the  reduced  mechanism. 


3.4.2  Variation  of  Burning  Velocity  as  a  Function  of  Pressure 

Calculations  were  done  to  study  flame  structures  with  changes  of  pressure  by  using  the  reduced 
mechanism  at  equivalence  ratio  being  unity.  The  burning  velocities  as  a  function  of  pressure 
is  shown  in  Fig.  3.26.  It  can  be  seen  that  the  burning  velocity  decreases  with  the  increasing  of 
pressure.  Converged  solutions  were  not  obtained  for  values  of  pressure  greater  than  1.4  atm. 
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Figure  3.20:  Variation  of  burning  velocities  of  premixed  CH4/NO2/O2  flames  calculated  using 
the  full  mechanism  and  the  reduced  mechanism  at  p  =  1  atm. 
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Figure  3.21:  Profiles  of  T,  H2O,  NO,  NO2,  CH4,  and  O2  calculated  using  the  fuU  mechanism 
and  the  reduced  mechanism  at  and  p  =  \  atm. 
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Figure  3.26:  Variation  of  burning  velocities  of  premixed  CH4/NO2/O2  flames  as  a  function  of 
pressure  calculated  using  the  fuU  mechanism  at  <?!>  =  1 

3.5  The  Reduced  Mechanism  for  Premixed 
CH2O/NO2/O2  Flames 

All  the  calculations  in  this  section  are  done  for  flame  systems  with  the  first  definition  of 
equivalence  ratio,  (/>,  and  at  pressure,  p  =  1,0  atm.  Two  cases  are  considered.  For  the  first 
case,  the  steady-state  assumptions  are  made  for  four  species  H2O2,  HCO,  N  and  HO2.  For 
the  second  case,  six  species  H2O2,  HCO,  N,  HO25  0?  OH  are  assumed  in  steady-state. 
The  burning  velocities  as  functions  of  equivalence  ratio  (f>  are  plotted  in  Fig.  3.27.  The 
results  corresponding  to  the  first  case  agree  with  those  calculated  using  the  fuU  mechanism. 
However,  the  results  corresponding  to  the  second  case  deviate  from  those  calculated  using  the 
fuU  mechanism,  especiaUy  in  the  vicinity  of  stoichiometry. 

3.6  Conclusions 

Numerical  calculations  were  performed  to  study  the  flame  structures  of  premixed  CH4/NO2/O2 
and  CH2O/NO2/O2  flames.  Two  kinds  of  definition  of  equivalence  ratio  were  defined  to  study 
the  influences  of  equivalence  ratio  on  burning  velocities.  The  flame  systems  with  the  first  kind 
definition  behave  like  a  usual  two-reactant  (fuel  and  oxidizer)  flame  systems.  The  structures 
of  the  flames  were  plotted  to  study  the  combustion  processes  in  the  flame.  One  outstanding 
feature  of  CH4/NO2/O2  and  CH2O/NO2/O2  combustion  systems  is  that  NO2  acts  like  a  fuel, 
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and  competes  for  radicals  with  fuel  in  the  combustion  processes  and  thus  decreases  the  con¬ 
centrations  of  the  radicals j  which  makes  the  flames  weaker.  For  a  pure  CH4/NO2  premixed 
flame,  although  numerically,  it  is  possible  to  get  a  stable  flame,  this  flame  is  very  weak.  For 
a  pure  CH2O/NO2  premixed  flames,  the  burning  velocity  is  very  small  compared  with  those 
near  stoichiometric  conditions,  but  it  is  still  of  the  same  order  of  the  burning  velocities  of 
methane-air  flame  near  stoichiometry.  This  explains  the  difiiculties  encountered  in  stabilizing 
pure  CH4/NO2  premixed  flames  in  the  laboratory. 

The  structures  of  CH4/NO2/O2  and  CH2O/NO2/O2  premixed  flame  systems  resembles 
those  of  hydrocarbon  flames.  They  include  a  preheat  zone,  an  inner  layer  and  a  post-flame 
zone.  The  inner  layer  can  further  be  subdivided  into  a  fuel  and  oxidizer  consumption  layer  in 
which  the  fuel  and  oxidizer  are  consumed  to  form  CO,  H2,  NO  and  some  H2O,  an  oxidation 
layer  in  which  CO,  H2,  and  NO  are  oxidized  to  form  CO2,  H2O  and  N2. 

A  seven-step  reduced  mechanism  is  developed  for  CH4/NO2/O2  premixed  flames.  The 
burning  velocities  calculated  using  the  reduced  mechanism  are  in  agreement  with  those  cal¬ 
culated  from  the  full  mechanism.  The  temperature  profile,  the  profiles  for  species  H2O,  NO2, 
CH4,  O2,  NO,  CO,  CO2,  H2,  N2  and  OH,  H,  0  predicted  by  the  reduced  mechanism  agree 
with  those  obtained  using  the  full  mechanism.  For  species  H2O2,  HCO,  HO2,  N,  CH3,  CH2O, 
and  CH3,  there  are  differences  between  the  values  obtained  using  the  reduced  mechanism  and 
those  from  the  full  mechanism. 

Two  sets  of  reduced  mechanisms  for  CH2O/NO2/O2  premixed  flames  are  tested.  The  first 
reduced  mechanism  is  obtained  by  assuming  four  species  HCO,  H2O2,  N  and  HO2  in  steady- 
state,  the  second  reduced  mechanism  is  obtained  by  assuming  six  species  HCO,  H2O2,  N,  HO2, 
0  and  OH  in  steady-state.  For  the  first  reduced  mechanism,  the  burning  velocities  predicted 
using  the  reduced  mechanism  agree  well  with  those  calculated  from  the  full  mechanism,  and 
the  temperature  and  species  profiles  for  the  stable  species  and  radicals  obtained  both  from 
the  reduced  and  the  fuU  mechanism  are  in  good  agreement.  For  intermediate  species  HCO, 
HONO,  HNO  and  HO2,  the  results  from  the  fuU  mechanism  and  the  reduced  mechanism  are 
also  in  agreement  with  each  other.  For  species  H2O2,  the  peak  value  obtained  from  the  re¬ 
duced  mechanism  is  of  two  orders  larger  than  that  calculated  from  the  full  mechanism. 
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Figure  3.27:  Burning  velocities  of  premixed  CH2O/NO2/O2  flames  calculated  using  the  full 
and  the  reduced  mechanisms  at  p  =  1  atm. 
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